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Abstract 



We consider the linear growth-fragmentation equation arising in the modelling of cell division 
or polymerisation processes. For constant coefficients, we prove that the dynamics converges to 
r"| . the steady state with an exponential rate. The control on the initial data uses an elaborate L 1 - 

norm that seems to be necessary. It also reflects the main idea of the proof which is to use an 
anti-derivative of the solution. The main technical difficulty is related to the entropy dissipation 
rate which is too weak to produce a Poincarc inequality. 
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■ 1 Introduction 

O 

We consider the growth-fragmentation (or cell division) equation 

% . { dt - ' dx ' J ™ - 

n(x = 0,t) = 0, 



d d f°° 

— n(x, t) + Q^n(x, t) + kBn(x, t) = B J k(x, y)n(y, t)dy, 



together with an initial condition n(x,t = 0) = n°(x) € L 1 (0,oo). Here, the number k > 1 represents 
the average number of pieces after fragmentation, while B > measures the relative intensity of 
the fragmentation process versus the growth process and x i— > k(x, y) gives the size repartition after 
fragmentation of a polymer (or division of a cell) of size y. Throughout this paper, the fragment 
distribution function k is assumed to enjoy the following properties: 

k(x, y) = for < y < x, (2) 

y ry 
K(x,y)dx = k, / xn(x,y)dx = y. (3) 
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These assumptions are usual in modelling fragmentation phenomena [H] and lead to two balance laws 
for the number of fragments J °° n(x, t)dx and the total mass Jq°° x n(x, t)dx, namely 

poo roo 

/ n(x,t)dx = Q := / n°(x)dx for all t > 0, (4) 
Jo Jo 

j roo roo roo 

— / x n(x,t)dx + (k — 1)B / xn(x,t)dx= / n(x,t)dx = g for all t > 0. (5) 
«* Jo 7o Jo 

These are instrumental for the existence theory of weak solutions. 

The evolution equation (pQ) is already renormalised so that zero is the first eigenvalue of the under- 
lying linear integro-differential operator. Indeed, recall that, under fairly general assumptions on re 
that were stated for instance in [HUH], there exists a unique steady state TV" € L 1 (0, oo) satisfying 

-k-N{x) + kBN{x) = B / k{x, y)N(y)dy, 

roo ( 6 ) 

N(x = 0) = 0, N > 0, / N(x)dx = 1. 

Jo 

In addition, convergence of solutions n to ([1]) with initial number of fragments g towards gN is estab- 
lished in \12\ Section 4] under suitable assumptions on k but without providing a rate of convergence 
(see also \14:\ Chapter 4]). These convergence results are based on the analysis of the entropy dis- 
sipation that we recall in Section [2j However, an estimate on the time needed to reach the steady 
state is of importance: in particular, an exponential rate of convergence guarantees that the steady 
state is reached fast enough to be observed in practice. Another motivation to investigate the rate of 
convergence is the study of the stability of steady states in nonlinear cases, at least in weakly nonlinear 
regimes. Such nonlinear models are met for example in prion dynamics [21 [9], cell division with 
finite resources [10} 114]. network communications [T], and many other fields. 

There are two particular cases where an exponential decay is proved. The first is equal mitosis 
equations (each cell divides in two cells of the same size) which corresponds to the choice k(x, y) = 
A5(y — 2x) (see Section [2]), and the method of proof in [HI [15] is extended here. The second is the case 
of age structured equations which is very particular and more classical because it corresponds to usual 
modelling; it corresponds to the formal limit a — ► in the case of general mitosis below in Section [2] 
and the exponential decay rate is known even the case of non-constant coefficients, see [4l \7\ llO]. 

The aim of this note is to show that such an exponential decay is actually valid for a broader class 
of fragment distribution functions re, namely those satisfying the positivity condition 

0(x,y):=~J K(z,y)dz>0, (7) 

in addition to ([2]) and ([3|). Below we give several examples for which these properties are met. The 
interest of our result is not only to extend the exponential trend to equilibrium to a wider class of 
data but also to explain the general structure by the positivity property in ([7]). The semi-norm on the 
initial data that comes out from our analysis is defined by 

roo fiyr roo poo 

HHII := / —{x) + kBM(x)-B /3(x,y)M(y)dy dx + (k + l)B \M(x)\dx (8) 
Jo dx Jo Jo 

with M(x) = Jq(ti(z) — g N(z))dz and g = J °° n(y)dy (following the notation introduced in dU). 
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Theorem 1.1 With assumptions ([3p, ([?]) and an initial condition satisfying \\\n \\\ < oo, the 
corresponding solution n to (CP satisfies 

\n(x,t) - gN(x)\dx <\\\n°\\\e- {k ' 1]Bt for all t>0. (9) 



Up to our knowledge it is not possible to obtain exponential decay with only control of L p -norms on 
the initial data (even though we are not aware of an explicit counterexample). This is related to the 
main difficulty in analysing the decay rate; the entropy dissipation is too weak to produce a Poincare 
inequality. We explain this in more details in Section [2j This is usual and the case of Fokker-Planck 
equations and "hypocoercive" equations are other examples [3| [16] . 

The outline of the paper is as follows. In the next section, we give some examples where the 
assumptions (J2]), ©, ([7]) are satisfied. Then we prove Theorem ll.il The proof uses an integral change 
of unknown which has already been used in [U [T¥| [T5] for similar purposes (see also [5j El [13] where 
it is applied to the uniqueness issue for related models). More precisely, we study the equation solved 
by an anti-derivative of the solution to ([1]). A crucial point here is that it satisfies a closed equation, 
a property due to the constant fragmentation coefficient kB. We gather several remarks in the last 
section. 

2 Examples and entropy dissipation 

We present here several standard examples of coefficients k(x, y) that our method is able to handle. 
We also recall the relative entropy associated with the equations at hand and indicate why a Poincare 
inequality cannot hold true. 

Equal mitosis. The most general model describing cell division relies on the choice 

k(x, y) = A5{y - 2x) = 25{x - y/2), k = 2, (10) 

that is, each cell divides in two cells of equal sizes. Assumptions ([2]) and (J3]) are obviously satisfied. 
Furthermore, we can compute 

X 

k(z, y)dz = 21 {x > y/2 } = 21 {y < 2x }, y) = 25(y - 2x) = 5(x - y/2), 

and (J7]) is satisfied. 

Uniform fragmentation. This is the simplest and most standard model for (physical or biological) 
polymers. It assumes uniform distribution of the fragments after each fragmentation event, i.e., 

2 

K(x,y) = -li x<y \, k = 2. (11) 

y 

Assumptions ([2]) and ([3]) are again obviously satisfied. In addition, 

[ n(z,y)dz = 2 mm { x >y} ^ 0(x,y) = 2 >x} , 

Jo y y 
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and is also satisfied. 



General Mitosis. An extension of equal mitosis is to take, for some parameter a with < a < 1/2, 

k(x, y) = -S(y--) + -^—5 (y - — ^— j = 5(x - ay) + 5(x - (1 - a)y). (12) 
a \ a J 1 — cry 1— cry 

Assumptions (|2|) and ([3]) are clearly fulfilled, still with k = 2, and 

K(z,y)dz = l{ x >ay} + ^{x>(l-a)y} = ^{y<x/a} + %{y<x/(l-a)} j 

y) = 5 - ^ + (5 h/ - 73-^ = aS ( x ~ a v) + i 1 ~ a ) S ( x ~ i 1 ~ a )y)- 
Therefore (|7|) is again satisfied. 

Homogeneous fragmentation. This generalisation of uniform fragmentation consists in using 

x a 2 + a 

K(x,y) = (2 + a) — l^j, k = YT^- (13) 

for some a > — 1 and satisfies fl2J) and ([3]). Also 

min{x,?/} 1+a i 1 " 1 "" 
k(z, y)dz = k , (3{x, y) = (2 + a )^_i {l/>a . } , 

and is again fulfilled. 

All these examples illustrate the main difficulty in order to prove the exponential decay of solutions 
and why entropy methods are not enough. We recall from |12^ [T4"l [16] and the references therein that 
the relative entropy inequality for equation (JTJ) reads 

d r * T , ^f n M 



dt 



L°° n{x)h (w) dx =- D ^ n w 



dx dy. 



,oo r / / n \ rT /n(x)\ rW / n(x) \ / n(y) nix) 

For a convex function H this is a non-negative quantity and this entropy inequality is usually enough 
to prove the strong convergence of n(t) to gN in a norm that depends on the choice of H. However 
the simplest method to prove exponential decay is through a Poincare inequality (whenever available) 
which asserts the existence of v > such that 

for all functions n satisfying n(x)dx = 0. 
The usual setting for the quadratic entropy H(r) = r 2 is to say that u = n/N satisfies 

/■oo poo roo 

v I N(x)u(x) 2 dx < D2[u] := B / / n(x,y)N(y)(u(x) — u(y)) dx dy, 
Jo Jo Jo 
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whenever J °° u(x)N(x)dx = 0. As far as we know, the validity of such an inequality is unclear and it 
is actually false in the case of equal mitosis. Indeed, in that case, k is given by (jlOp and the right-hand 
side D2M is 

/•oo 

D 2 [u] = AB N(2x)(u{x) - u(2x)) 2 dx, 
Jo 



which clearly vanishes for the class of functions that satisfy u(x) = u(2x) for all x € (0, 00): such 
functions are obtained from a given bounded function (p defined on [1,2] by setting u{x) = 4>(x/2 :) ) 
for x G [2i,2i +l ) and j £ Z. 



3 Proof of Theorem 11.11 

As already mentioned, the proof uses an integral change of unknown which is our starting point. 
Step 1. (change of unknown function) Let n be a solution to ([1]) and set 

/•X 

M(x,t) = / [n(z,t) - gN(z)]dz for (x, t) £ [0, oo) 2 . 
Jo 



It satisfies M(0,t) = M(oo,t) = thanks to the conservation of the number of fragments (J2|) which 
implies 

/'OO poo roo 

[n(z, t) - gN(z)]dz = / n°(z)dz -pi N(z)dz = 0. 
Jo Jo 



The motivation for using M is that it satisfies a closed equation of the same type, namely, 



-|m(x, t) + ^M{x, t) + kBM(x, t) = B J 0(x, y)M (y, t)dy, 
M(x = 0,t) = 0. 



(14) 



Indeed, integrating the first equation in ([T]) over (0,x), x £ (0, 00), we obtain 



0_ 

Of 



d r f°° 

M(x,t) + —M(x,t) + kBM(x,t) = B / K,(z,y)n(y,t)dydz 

U x Jz=0 Jy=0 



y= 

B I j-M(y,t) / K(z,y)dzdy 
y =o oy J z=0 



B 



d_ 

dy 



k(z, y)dz 



z=0 



ly=0 

/•oo 

B / (3(x,y)M(y,t)dy. 
Jo 



M(y,t)dy 



Step 2. (Property of (5) To proceed further, we need a specific property of j3, namely 



f3(x,y)dx 



L^y{L K{z ' y)dz ) dx 



d_ 

dy 

d ( r y 



x=0 Jz=0 
V 



n(z,y)dzdx I + / n(z,y)dz 
J Jz=0 

v \ 
dxdz ] + k 



dy 
d_ 

dy 



z=0 



(y - z)n{z,y)dz + k 



and thus 



(k - l)y] + k 



P(x, y)dx = 1. 



(15) 



Step 3. (Exponential decay of M) We now argue on the function M. Owing to the positivity ([7]) of /3, 
we infer from (|14|) that 

-\M(x,t)\ + —\M(x,t)\+kB\M(x,t)\< / P(x,y)\M(y,t)\dy. 

at dx j 

After integration of the above inequality with respect to x over (0, oo), we deduce from (fT5j) that 

j rco poo poo poo 

-J \M(x,t)\dx + kB I \M(x,t)\ < I J (3(x,y)\M(y,t)\dydx 

poo 

< B I \M(x,t)\dx, 
Jo 



\M(x,t)\dx < e - {k - 1)Bt / \M(x,0)\dx for t > 0. 



and thus 



Step 4- (Exponential decay of n) To transfer the exponential decay to the solution n itself is not always 
possible, see for instance the coagulation case |S]. Here, we can follow [15] and argue as follows. We 
first notice that dM/dt satisfies the same equation as M. Therefore Step 3 and (Q3|) also give 



f 



dM 



dt 



(x,t) 



dx < e~^ Bt 



dM 



< e -{k-l)Bt 

Next, owing to ([7]) and ([To]) , we have 



dt 

dM 



(x,0) 



dx 



dx 



poo 

(x, 0) + kBM(x, 0) - B / P(x,y)M(y,0)dy 

Jo 



dx. 



P(x,y)M(y,t)dy 
As a consequence, (fT4j) also reads 

n(x, t) — gN(x) = - 



dx < 



\M(y,t)\ / P{x,y)dxdy 



\M(y,t)\dy. (16) 



dM 
~dt 



(x, t) - kBM(x, t) + / p(x, y)M(y, t)dy, 



6 



and we conclude with the help of (|16[) that 

dM 



oo poo 

\n(x,t) — gN(x)\dx < / 

Jo 



Of 



(x,t) dx + (k + l)B / |M(x,t)|dz 

JO 

/•oo .g/l/f /*oo 

< e-W* J —(x,0) + kBM(x,0)-Bj^ (3(x, y)M(y, 0)dy 

pOO 

+ {k+l)Be~ {k ^ Bt / |M(x,0)|di 



dr 



as claimed. □ 



4 Concluding remarks 

(a) Owing to (fT6j) and its definition in ([8]), the semi- norm actually satisfies 

dM 



n u < 



dx 



(x,0) 



roc 

dx + 2(k + l)B / \M(x,0)\dx, 
Jo 



and |||n°||| is thus dominated by the Vl^ 1 ' 1 -norm of M(.,0). 
(b) Theorem 11.11 can be extended to growth- fragmentation equations with non-constant velocity 

( d d f°° 

— n(x, t) + — (r(x)n(x, t)) + kBn(x, t) = B J k(x, y)n(y, t)dy, 

n{x = 0,i) = 0. 

The fragment distribution function k is still assumed to fulfil the properties ([2]), (J3|), ([7]), and 
the velocity r is required to satisfy 

dr 

T(x)>T m >0 d:={k-l)B- sup — (<e) > 0. 



xG(0,oo) 



9x 



Under these assumptions, the exponential decay rate is e Indeed, the equation on M is 

— M(x, t) + t(x)—M(x, t) + kBM(x, t) = B (3(x, y)M(y, t)dy, 
< at dx J 

M(x = Q,t) = 0. 

Therefore, proceeding as in the third step of the proof of Theorem II. II we have 

— I \M(x,t)\dx + (k-l)B \M(x,t)\dx= \M(x,t)\ — (x)dx. 
dt J Jo Jo ox 

The decay rate for M follows directly from this equality. The positive bound from below required 
on r is next used to transfer this decay rate on n, arguing as in the fourth step of the proof of 
Theorem II. II 
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(c) The exponential decay rate in Theorem 11.11 seems to be optimal: indeed, consider the case where 
k is given by (flOj) (equal mitosis), that is, k(x, y) = A5(y — 2x) = 25(x — y/2), and introduce the 
function £ defined by 

. . . IdN /x\ /n . 



2 dx \2 

Owing to the properties of TV given in © and the specific choice of k, it is straightforward to 
check that (x, t) i — > e~ Bt ^(x) is the solution to {!]) with initial condition £ and the integral of 
£ vanishes. Consequently, given g > 0, n(x,t) = gN(x) + e~ Bt ^(x) is the solution to (TTJ with 
initial condition n° = gN + £ and 

/■oo 

\n(x,t) - gN(x)\dx = e~ Bt / |£(x)|da; 

J o 

(recall that k = 2 in that case). 

(d) Similarly, when n{x,y) = 21^ x<y y/y (which corresponds to uniform fragmentation, see (fTTjl ). 
the equilibrium is known explicitly and given by N(x) = 4B 2 xe~ 2Bx for x > 0. Introducing 
?/(x) := x(fix — 2)e~ Bx for x > 0, one readily checks that (x,t) i — > e~ Bt r](x) is the solution to 
([T|) with initial condition r]. Observe that the integral of rj vanishes. As in the previous case, 
given q > 0, n(x, t) = gN(x) + e~ Bt r]{x) is the solution to (H|) with initial condition n° = gN + 77. 
Since 

roc poo 

\n(x,t) - gN(x)\ dx = e~ Bt / \rj(x)\dx, 



and k = 2, this example also yields the optimality of the decay rate obtained in Theorem II. II for 
that particular case. 
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